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Abstract. A coarse-graining strategy for dilute and semi-dilute solutions of 
interacting polymers, and of colloid polymer mixtures is briefly described. Monomer 
degrees of freedom are traced out to derive an effective, state dependent pair potential 
between the polymer centres of mass. The cross-over between good and poor solvent 
conditions is discussed within a scaling analysis. The method is extended to block 
copolymers represented as "necklaces" of soft "blobs" , and its success is illustrated here 
in the case of a symmetric diblock copolymer which exhibits microphase separation. 
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1. Introduction: a coarse-graining strategy 

The present paper summarises a collective effort by the Cambridge group and others 
to bridge the "cultural divide" between statistical mechanics of polymer solutions and 
melts on one hand and of "simple liquids" on the other. While the former is dominated 
by field-theoretic methods and scaling concepts pioneered by S.F. Edwards and P.G. de 
Gennes, the description of simple liquids, which lack the scale invariance of polymers, 
requires a more atomistic approach. To bridge the gap we have developed a systematic 
coarse-graining strategy first suggested by Flory and KrigbaumpP, whereby individual 
monomer degrees of freedom are traced out for fixed centre-of-mass (CM) positions of 
the polymer coils, thus defining an effective interaction between the CMs. Consider a 
system of N identical polymer, each consisting of M monomers (or segments) at positions 
'^ia i ^ N] 1 < a < M). If H{{ria}) is the interaction hamiltonian of the system, 
the probability distribution of the N CMs at positions Ri is 

1 „ N 1 N M 

PN{{m = ^ / e-/^^(^-» n m - tj e ^^^) n n (i) 

-/3v;;/({flj) 

= = ^ (2) 

where /3 = Qa? is the partition function for the N x M monomers, and the total 
effective interaction energy of the N CMs is rigorously defined by: 

K//({^.}) = -kBT\n[C X P^({i?.})] (3) 

This effective energy is a free energy, and hence state-dependent, and, in general, many- 
body in nature. In the low concentration limit, reduces to a sum of effective pair 
interactions between two isolated polymer coils: 

1^2(1^1 - R2\) = -ksTlnlC X P2{Ri,R2)] (4) 
Since global properties of polymeric systems are independent of chemical detail in the 
scaling {L —>■ oo) limit, we adopt henceforth a simple lattice model of polymer solutions 
namely that of mutually and self-avoiding walks (SAW) of length L = M — 1 on a 
cubic lattice of lattice spacing b (equal to the segment length); non-connected nearest 
neighbour monomers of the same or different polymer coils have an attractive energy 
— e. This model accounts for the key polymer features, namely connectivity, excluded 
volume and solvent quality (through the value of e). 

Convenient thermodynamic variables are the polymer density p = N/ {Qb^) {Q being 
the size of the lattice) , the monomer density c = Mp (equivalently the monomer packing 
fraction = cb^ equal to the number of lattice sites occupied by monomers) and the 
temperature T (/3* = e/ZcsT). A key characteristic is the overlap density p* = S/AnRg^ 
(where Rg is the radius of gyration, ~ bL'^) which corresponds to the cross-over from 
the dilute {p < p*) to the semi-dilute regimes (p > p*). The semi-dilute regime differs 
from the melt in that the monomer packing fraction remains negligible; for any given 
p/p* this is only achieved for sufficiently long polymers, since = (p/ p*)L^^~^'^^ ~ L~^^^ 
for SAW polymers. 
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2. Polymers in good solvent 

Consider first tlie athermal limit of SAW polymers (e = 0). In the low (p — > 0) 
density limit P2{Ri,R2) is then simply the probability that there is no monomer- 
monomer overlap for a fixed distance r = — i?2| between the CMs. This is well 
adapted to Monte-Carlo (MC) sampling. Scaling theory predicts that the resulting 
effective interaction at full overlap, V2{r = 0) is independent of chain length L^]. Early 
simulations with rather short chains pointed to V2{r = 0) ~ 2/cBTj3j. An in-depth 
investigation of the L-dependence shows that V2{r) is well represented by a Gaussian 
V2{r) ~ uexp{—a{r / RgY) , where a is of order 1, while u/ksT ^ 1.80 
At finite polymer concentration p, three and more-body effective interactions come into 
playjH]. A more efficient strategy is to determine an effective density-dependent pair 
potential f2(r;p) by inverting the CM-CM pair distribution g{r) from full monomer- 
level MC simulationsjH] [7|. It was proven that this inverse problem has a unique 
solution[Hj. The inversion procedure is implemented using the HNC-integral equation|nj. 
The Ornstein-Zernike relation [Hj allows the extraction of the direct correlation function 
c(r) from the MC data for h{r) = g{r) — 1 at any given density. The HNC closure then 
expresses f2(r) as [H]: 

[3v2{r) = — ln[(7(r)] + h{r) — c(r). (5) 

The first term on the r.h.s. is the potential of mean force, while h{r) — c(r) describes 
the effect of correlations. 

The MC generated g{r) show that correlations decrease as p increases, contrary to 
the more familiar behaviour observed for hard core systems. The overlap value g{r = 0) 
increases steadily toward 1 (the ideal gas value) confirming that in the high density limit 
of a melt, polymer chains indeed behave as non-interacting polvmers|lUj. The resulting 
effective pair potential is only moderately density dependent and is well fitted by a sum 
of gaussians. The range of V2{r) tends to increase with p, and the potential develops a 
small amplitude negative tail for r significantly larger than B.^ |^ [TT] . 
The link with thermodynamics is via the compressibility relation jH], which allows the 
osmotic pressure P to be expressed as: 

f3P{p)= r[l-p'c{k = 0;p')]dp' (6) 
Jo 

where c{k) is the Fourier transform of c(r). Use of the effective pair potential in 
conjunction with the virial and energy equations is meaningless jl2j . The equations of 
state calculated from MC simulations of the full monomer level polymer representation 
and from simulations based on the much less CPU-intensive effective potential 
representation agree within numerical uncertainties, underlining the adequacy if the 
HNC inversion procedure for such "soft" effective particles. Well into the semi-dilute 
regime (p ^ p*) the slopes of the calculated equation of state agrees with the des 
Cloizeaux scaling prediction PP ^ p^'^/i^'^-^) ^ p^^^, where u = 0.588 ~ 3/5 is the Flory 
exponent for the radius of gyration in good solvent ~ &L'^)PT]]. 



Polymer solutions 



4 



Neglecting the density- dependence of V2{r), i.e. extending the low density gaussian form 
to all densities, brings us back to the "Gaussian core model" (GCM) first introduced by 
Stillinger^lBj. which exhibits interesting behaviour at low temperatures (/?* 3> l)|14j. 
In the regime relevant for polymer solutions (/?* ~ 1) the model leads to "mean field 
fluid" behaviour at sufficiently high density, where the random phase approximation, 
c(r) = —[3v2{r) leads to a quadratic equation-of-state for p ^ p*^^^Ei- 

f3P = p+^(3v2{k = 0)p' (7) 

Incorporating the p-dependence of f 2 is thus seen to change the asymptotic behaviour 
into des Cloizeaux scaling p^^^. As a by-product of the GCM, we have developed 
a multiple occupancy lattice model, which also gives rise to interesting microphase 
separation [in]. 



3. From good to poor solvent conditions 

We now turn our attention to the case where adjacent monomers attract, i.e. e 7^ or 
/?* > 0. This attraction is solvent induced, and the quality of the solvent deteriorates 
as P* increases, leading to contraction of the polymer coils. At the 9 temperature {Pq = 
e/kBTg), repulsion and attraction between polymers cancel, at least in the low density 
limit, so that polymers exhibit the scaling behaviour of ideal polymers(/2g ~ L^/^). 
Below Tg, polymer coils collapse into globules (Rg ~ L^^^), and phase separation occurs 
into polymer-rich and polymer-poor solutions. 

The most convenient diagnostic for locating Tg from simulations is to calculate the 
second virial coefficient B2{L] T) as a function of temperature and polymer length. The 
L-dependent Boyle temperature Tb{L) is that at which B2{L] T) vanishes for a fixed L, 
then: 

Tg = hm Tb{L) (8) 

L— +00 

This leads to the estimate = 0.2690 ±0.0002 Note that B2{L; T) can be directly 
expressed in terms of the low density limit of the effective CM pair potential: 

B2{L; T) = 27r / [1 - e-^^^^^'^'^V'^^^ (9) 
Jo 

Extensive MC simulations were used to determine the effective pair potential for fixed 
length L = 100 over a wide range of temperatures(0 < f3* < 0.3) and densities 
p[TS]. As P* increases, V2{r = 0) decreases and V2{r) develops an attractive tail for 
r > Rg. Eventually V2{r) violates Ruelle's necessary condition for the existence of a 
thermodynamic limit, namelv[T^: 

l2 = J V2{r)df> (10) 

Since V2 depends on L, p and T, so does I2, and for any given p and L, the limit of 
stability temperature Tg is determined by 

/2(L,p,r = r,) = (11) 
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Below Tg, single polymer coils will collapse and we conjecture that Tq = 
lim — .n lim t ^^ TJ p. D [TS] . However, if p is increased at a given temperature T < Tg 
the effective pair potential increasingly reverts to being repulsive until the Ruelle 
criterion PJ is satisfied, and the polymer solution becomes thermodynamically stable 
again. This "restabilisation" reflects a phase separation scenario under poor solvent 
conditions(T < Tg). 

The variation of V2{r) with temperature and density, as extracted from the HNC 
inversion of MC data, is semi-quantitatively reproduced by solutions of the PRISM 
integral equation|2ni for the thread model(6 — 0, L ^ oo at fixed Rg) of polymer 
solutions[2I]|22]- PRISM assumes all monomers to be equivalent {i.e. neglects ends 
effects) and yields the monomer pair distribution function gmm{r) of a polymer solution. 
The CM distribution function g{r) required to extract the effective pair potential V2{r) 
may be related to gmm{r) by an approximate, but accurate relation involving internal 
form factors of a single coil [231 • Further progress along these lines might eventually 
bypass the need for time-consuming simulations of full monomer level models required 
to determine g{r). 

The equation of state may be computed as a function of p/p* and /?*, using a somewhat 
cumbersome method based on the contact theorem |24j |25j. or much more efficiently by 
subjecting the polymers to a gravitational field and invoking hydrostatic equilibrium |26j. 
If p{z) denotes the CM(or monomer) density profile of the polymers in a vertical field, 
which is easily measured in simulations, the osmotic pressure at an altitude z is simply 



where ( = ksT/Mg is the gravitational length, which must be chosen (by tuning the 
product Mg) significantly larger than Rg for the macroscopic description (which follows 
from the local density approximation within density functional theory of inhomogeneous 
fluids) to hold[2Zj. Elimination of the altitude z between P{z) and p{z) then leads to 
the bulk equation of state P{p). Examples from simulations of L=500 chains are shown, 
for four temperatures in figure 1. 

4. Long polymers: corrections to scaling 

Polymers are critical objects in the universality class of the zero component (n=0) limit 
of the n- vector model of critical phenomena |28j. The L — oo limit is equivalent to 
the limit of divergent correlation length ^ at the critical point of a second-order phase 
transition. The properties of very long chains can hence be investigated using the 
powerful method of the renormalization group (RG) and field theory |29|. These predict 
two scaling regimes, an athermal one corresponding to the good solvent (SAW) limit, 
and the second corresponding to the ^-solvent regime where ideal polymer statistics hold. 
The cross-over between the two regimes is discontinuous in the scaling limit L — oo. 
The objective is to predict the behaviour for large but finite L from a finite size scaling 
analysis of MC data[3n||S]- We restrict the discussion to the good solvent regime. For 




(12) 
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Figure 1. Equation of state on a log-log scale, calculated using Dickman's (contact 
theorem) method, (symbols) and hydrostatic equilibrium(solid lines). 



l3 < Pe any universal (dimensionless) ratio R may be represented as 

R{L,P) = R* + ^^ + ... (13) 

where R* is the temperature independent scahng hmit [L oo) and the exponent 
A = 0.517jni]. Higher order terms involve exponents of the order of 1 or larger. The 
temperature-dependent coefficient aR{(3) is non-universal, but ratios of such coefficients 
for two different dimensionless qualities are again universal, i.e. model-independent. 
A much studied example of a universal ratio is A2{L,[3) = B2{L, P)/ R^{L, P). A 
similar ratio involving the third virial coefficient A^lL, /5) = B^{L, /5) / Rg{L, /5) has been 
examined in detail in reference 0. is found to be always positive, but to go through 
a sharp minimum near T = Tg. The effective pair potential between the CMs of two 
isolated polymers (p <^ p* limit), divided by A^^T, is another universal ratio, which has 
been investigated by a similar scaling analysis in [S]: 

Pv2{r,L,P) = v^{x) + ^v,{x) + ... (14) 

where x = r/Rg. v^o, fc and a^lP) are extracted from a careful analysis of MC data 
for several lengths L and inverse temperatures < P < Pg. These may then be used to 
predict Pv2 for any length and temperature, and agreement with MC data is excellent 
for P* < 0.2 and L > 500. Closer to the ^-temperature the convergence of ()14|) with L 
is found, not surprisingly, to be much slower, and one must then switch to the scaling 
analysis appropriate for the 6'-regime. The scaling analysis has been recently extended 
to binary mixtures of polymers of different lengths and to star polymers [22] 



5. Soft polymers and hard colloids 

Mixtures of colloidal particles and non-adsorbing polymers have attracted considerable 
experimental and theoretical attention over the last two decades, because of interesting 
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phase behaviour induced by the famihar depletion mechanism |33j. The effective 
depletion interaction between spherical colloidal particles may be tuned by varying the 
size ration q = Rg/Rc (where Rc is the colloid radius), and the polymer concentration 
p. Consider first the case of polymers between 2 plates(g = 0) separated by z. The 
depletion potential per unit area W{z) is defined by the difference in polymer grand 
potentials: 

W{z) = ]-[Q{z) -Q{z = oo)] (15) 

POO 

= / [P{z') - P{z = oo)]dz' (16) 

Clearly since at contact 2 depletion zones are destroyed W{z = oo) = —2'yw{p), where 
7iy is the polymer -wall surface tension. The simple ansatz: 

Wiz) = WiO) + Pip)z ■,z<DM = -^^ (IT) 

= ;z>DM (18) 

reproduces direct simulation data well. The result for SAW polymers |S| differs 
considerably from that for ideal polvmers|33j: the range Dw is shorter for interacting 
polymers, and decreases with increasing density, while it is independent of density for 
ideal polymers; the contact value W{0) decreases faster with p for interacting polymers 
since the surface tension scales as p^^'^ in the semi-dilute regime while it scales as 
p for ideal polymers. For finite q < 1, the depletion force between two spheres can be 
approximately related to the depletion potential between 2 planes via the Derjaguin 
approximation, and by correcting for the decreasing range of the force due to partial 
wrapping of the polymer coils around the spherical colloid [Mj. 

Such pair depletion interactions do not, however, account for effective many-body 
interactions due to finite colloid concentrations. To that purpose the coarse-graining 
strategy described in sections 1 to 3 may be extended to the two component colloid- 
polymer system. An effective state-dependent colloid-polymer pair interaction Vcp{r) 
may be extracted by an HNC inversion of the polymer density profile around a 
sphere ^I] [SI], similar to the inversion procedure used to determine the effective 
polymer-polymer pair potential fpp(r). The colloid-colloid pair potential Vcc{r) is well 
approximated by a simple hard sphere interaction. MC simulations of this effective 
two-component system were used to calculate the phase diagram of the mixture for 
q < I [Sni- The calculated binodal agrees very well with experimental datajSH] for 
interacting polymers. Significant qualitative differences arise between phase diagrams 
for idealj23 ^ind interacting polymers, particularly for larger q: the range of the 
concentrated ("liquid") colloidal phase in the colloid density-polymer density plane is 
considerably reduced when polymer interactions are included, and the critical point 
occurs at significantly higher packing fraction of the two species|2n]- These trends 
become more pronounced in the "protein limit" of large polymers and small colloids 
(g 3> l)|Sn]- Solvent quality has a strong influence on the induced depletion interaction 
between colloids, with polymers under 6 conditions leading, not surprisingly to a pair 
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interaction close (but not identical) to that induced by ideal polymers |25j. at least for 
p/p* < 1. 

6. Diblock copolymers and beyond 

The coarse-graining strategy may be extended to polymers other than linear 
homopolymers considered so far. Star polymers, for instance, have been investigated 
along similar lines [10]; the CM is replaced by the midpoint where the / arms of a star 
polymer meet; the resulting effective pair interaction diverges logarithmically as r ^ 
and hardens as / increases. 

The case of symmetric diblock copolymers AB has been examined very recently jl^. 
The A and B strands are represented as soft "blobs", the CMs of which are tethered 
by an "entropic spring" characterised by the intramolecular potential 4>AB{f) which is 
derived from the MC-generated distribution function of relative distances of the A and B 
CMs on the same copolymer. There are now three intermolecular CM potentials VAA{f), 
vbb{t) VAB{f)i in addition to the intramolecular potential. The inversion procedure to 
go from the partial distribution functions Qapi.^) to the pair potentials Vo,fj{r) is now 
more involved |41j. Even in the low density limit one already faces a four-body problem. 
The inversion procedure has been carried out in that limit for a simple athermal 
model, the ISI model, where A-A and B-B pairs behave like ideal polymers, i.e. freely 
interpenetrate, while A-B pairs behave like mutually avoiding walks. This is the block 
copolymer equivalent of the familiar Widom-Rowlinson model [42] which drives phase 
separation of simple atomic fluids (where A and B are untethered). Because the two 
strands of the AB block copolymer system are tethered, macroscopic phase separation 
due to A-B incompatibility is suppressed and reduces to microphase separation: the 
symmetric block copolymers form a lammelar phase, which was indeed observed in MC 
simulations of both the full monomeric and the coarse-grained representations!^ . The 
resulting equation of state calculated by the hydrostatic equilibrium method is shown 
in Figure 2. Z = [3P/p is seen to first increase linearly up to p/p* ~ 2, where it 
flattens out, and thereafter decreases slowly to give an asymptotic value > 1. This 
may be understood by noting that in the lamellar phase(which develops for p/p* > 2), 
the repulsive A-B contacts are greatly reduced. Figure 2 also shows the equation of 
state calculated (with much less computational effort) for the coarse-grained "soft- 
dumbbell" model with effective intra and inter-molecular pair potentials (pABi^) and 
Vajsir) determined in the zero density limit. The agreement is excellent up to p ~ p*, 
and remains semi-quantitative thereafter, despite the fact that the density dependence 
of the effective potential has not been taken into account. 

The "soft dumbbell" representation of diblock copolymers provides a hint of how to 
extend the coarse-graining strategy of linear homo or hetropolymers over a wide range 
of polymer concentrations. As the ratio p/p* increases, the fundamental length scale 
gradually crosses over from Rg (for p/p* < 1) to the correlation length ^ ~ p^^/^ deep 
in the semi-dilute regime, to the segment length b in the melt. The coarse-graining 
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Figure 2. Equation of state generated from MC simulations for diblock copolymers 
using full monomer level and coarse-grained "blob" models. 

strategy put forward in this paper applies to dilute and initial semi-dilute regime, where 
polymer coils are well represented by a single "soft-core" particle with a radius of the 
order Rg. Deeper into the semi-dilute regime, the blob picture|lSI applies for polymers 
confined by other polymers or in a pore; each polymer reduces to a "necklace" of blobs 
of radius C,, tethered by entropic springs; different blobs on the same or neighbouring 
chains interact via a quasi-gaussian soft-core potential, as introduced in earlier sections. 
The radius of each blob decreases, and hence their number increases (for a given overall 
length L) as the ratio p/p* increases, until the melt regime is reached where the blob 
size reduces essentially to the segment length b, and the coarse-graining strategy is no 
longer of any use. Over the whole semi-dilute regime a polymer may thus be pictured 
as a necklace of blobs, and the present coarse-graining strategy allows, in principle for 
an unequivocal determination of intra and inter-molecular effective interactions between 
these blobs. Work along these lines is in progress. 
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